Exotic Charges, Multicomponent Dark Matter and Light Sterile Neutrinos 
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Generating small sterile neutrino masses via the same seesaw mechanism that suppresses active 
neutrino masses requires a specific structure in the neutral fermion mass matrix. We present a 
model where this structure is enforced by a new U(l)' gauge symmetry, spontaneously broken at 
the TeV scale. In order not to spoil the neutrino structure, the additional fermions necessary for 
anomaly cancellations need to carry exotic charges, and turn out to form multicomponent cold dark 
matter. The active-sterile mixing then connects the new particles and the Standard Model — opening 
a new portal in addition to the usual Higgs- and kinetic-mixing portals — which leads to dark matter 
annihilation almost exclusively into neutrinos. 

PACS numbers: 14.60.Pq, 14.60.St, 95.35. +d, 14.70.Pw 



INTRODUCTION 



Neutrino oscillation experiments have by now firmly established the existence of neutrino oscillations 
and lepton flavor mixing, indicating that the Standard Model (SM) framework in particle physics has to 
be extended to include neutrino masses. Apart from the traditional three neutrino oscillation picture, the 
LSND Q and MiniBooNE Q short-baseline experiments suggest the presence of sterile neutrinos at the 
cV scale, which do not participate in the weak interaction but mix with active neutrinos with a mixing 
' r- 1 ' angle 6 8 ~ O(0. 1) 0,01 • Furthermore, recent re-evaluations of reactor anti-neutrino fluxes indicate that 

Q_[ the previous reactor neutrino experiments had observed a flux deficit, which can in fact be interpreted 

by additional sterile neutrinos with masses at the eV scale (see Ref. Q for an exhaustive overview of the 
field). Moreover, the light-element abundances from precision cosmology and Big Bang nucleosynthesis 
favor extra radiation in the Universe, which could be interpreted with the help of one additional sterile 
neutrino, albeit with a mass below eV 0. 

From the theoretical point of view, the question then arises how the small sterile neutrino mass scale 
0(eV) can be motivated compared to the electroweak scale of 0(100 GeV). Recall that the seesaw 
qq mechanism is one of the most popular theoretical attempts to understand the smallness of active 

. neutrino masses and also the baryon number asymmetry of the Universe. An obvious ansatz is then 

I/"") ■ to use the same seesaw mechanism to suppress sterile neutrino masses. To this end, the right-handed 

neutrino content has to be extended compared to that in the simplest type-I seesaw mechanism, and a 
specific flavor structure, i.e., the minimal extended seesaw (MES), has to be employed [1,0. Explicitly, 
in the MES model, the SM fermion content is extended by adding three right-handed neutrinos vri (for 
O^l ■ i = 1,2,3) together with one singlet fermion S,^ while the full Majorana mass matrix for the neutral 

fermions in the basis (u e , v^, v T , v R 1; v R 2 , v R 3 , S c ) is given by 



m D 

A^mes = I m T D M R m s I ■ (1) 



X: 

c3" V m s ° 

After integrating out the heavy right-handed neutrinos vrj, one obtains three massive eigenstates with 
masses around Mr, three with masses around the eV scale and one massless neutrino. Taking a hier- 
archical spectrum mo < < Mr, one of the three massive eigenstates will be a sterile neutrino and 
its admixture with active neutrinos is suppressed by a factor 0(ran/ms) (see Ref. @ and Sec. IIII Al for 
detailed discussions). 

Note that the MES structure defined in Eq. (TT|) can be obtained with discrete flavor symmetries 0], 
under which the right-handed neutrinos and S carry different charges. In addition, the MES structure 
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t To clear up potential confusion right away: all the Ufjj and Sj introduced in this paper are just right-handed fermions, 
sometimes referred to as singlets. We denote the Sj with a different symbol than vjh to emphasize that they are not the 
usual right-handed neutrinos from the seesaw mechanism, because they carry additional (hidden) quantum numbers and 
do therefore not partner-up with the active neutrinos. 
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could also be obtained in models with abelian symmetries. For example, one may introduce an extra U (1)' 
symmetry, and take the three right-handed neutrinos vrj to be neutral under U(l)' (all SM particles are 
neutral too). One may then write down a bare Majorana mass matrix Mr for ur, which is unprotected 
by the electroweak or Z7(l)' scale. The right-handed singlet S on the other hand carries a U (1)' charge Y', 
and we further introduce an SM singlet scalar cj> with charge —Y'. The gauge invariant coupling S c Vr4> 
then generates the mj matrix in Eq. ([1} after <f> acquires a vacuum expectation value (VEV), while the 
Majorana mass for S (i.e., S C S) and a coupling to the active Vl are still forbidden by the U(l)' symmetry 
at the renormalizable level. Such a simple realization of MES however suffers from the problem of triangle 
anomalies, and can therefore only work as a global symmetry, whose spontaneous breaking would result 
in a masslcss Goldstone boson. While this might not be disastrous, more interesting phenomenology 
arises when the U(l)' is promoted to a local symmetry. Consequently, one has to extend the model by 
additional chiral fermions so as to cancel the arising gauge anomalies. Along these lines, possible model 
constructions for sterile neutrinos in the U(l)' framework have already been discussed in Refs. using 
an effective-field-theory approach. 

In this note we will work in the seesaw framework and discuss minimal renormalizable and anomaly- 
free U(l)' models, which are spontaneously broken by just one additional scalar and reproduce the MES 
structure accounting for the 3 + 1 or 3 + 2 scheme of light sterile neutrinos. In particular, we will show 
that the additional singlet fermions employed for the anomaly cancellation turn out to be stable — due 



to remaining Zjv symmetries 11] — implying that they can be viewed as good candidates for the dark 



matter (DM) in the Universe. The remaining parts of this work are organized as follows: In Sec. [TT| we 
explain the framework of exotic charges under a gauged U (1)' and why they are needed to obtain the MES 
structure. In Sec. IIIII we discuss in some detail the phenomenology of a specific example with one light 
sterile neutrino (3 + 1 scheme) and three stable DM candidates, with a focus on the novel effects inherent 
in our model. We briefly discuss other interesting examples of this framework in Sec. IIVI including an 
extension to the 3 + 2 case. Finally, we summarize our work in Sec. [V] 



II. DARK SYMMETRY 



As already mentioned in the introduction, adding just one extra right-handed singlet S to the three ur's 
results in triangle anomalies if only S is charged under the extra U(l)' symmetry. Instead of treating 
U(l)' as a global symmetry, we gauge the U(l)' in the rest of the work, and accordingly introduce 
additional singlet chiral fermions to cancel the anomalies. As we will see below, these new states need 
to decouple from the neutrino sector in order not to spoil the MES structure, and automatically lead to 
DM candidates without the need for additional stabilizing discrete symmetries. 

For a gauged U(l)' symmetry under which all SM particles are singlets, there are no mixed triangle 
anomalies, so anomaly freedom reduces to the equations 

]>>'(/) = and £(y'(/)) 3 =0, (2) 

/ / 

where / stands for our new right-handed fermions. In order to cancel the contribution from S = Si, more 
U(l)' charged chiral fermions Si>2 have to be introduced. The solutions of Eq. (|2|) for n = 2 are simply 
given by Y'(S±) = —Y'(S2). In this case, a bare mass term mSfS^ — unconstrained by any symmetry — 
can be constructed, which spoils the desired MES structure for light sterile neutrinos unless we make 
m small. There is no integer solution with F' /0 for n = 3 according to the famous Fcrmat theorem. 
In the case of n = 4, it is easy to prove that there is no phenomcnologically interesting solution since 
two of the Si must have U(l)' charges of opposite sign and equal magnitude, inducing an unconstrained 
bare mass term as in the case of n = 2. For n > 5 however, there exist interesting non-trivial anomaly- 
free charge assignments — dubbed exotic charges hereafter — for example the set (10, 4, —9, 2, —7) for 
n = 5 [lOl Il2j . In order to make all new fermions massive at tree level with just one scalar 0, even more 
chiral singlets have to be introduced. For the 3 + 1 scheme discussed in the main text, we add seven 
singlet fermions Si to the model (the 3 + 2 scheme discussed in Sec. IIV Al needs six). The charges of all 
the ten right-handed fermions are listed in Tab. HI they are of course by no means unique, but serve as 
a simple illustration of this framework. We further stress that at least three U(l)' singlet right-handed 
neutrinos vrj are needed in order to explain the observed light neutrino mass-squared differences Amfi, 
Amjj and Am^, resulting in one massless active neutrino. This is however not a hard prediction of the 
MES scheme; adding a fourth Vr (or even more) to the model makes all light neutrinos massive and does 
not qualitatively change or complicate the discussion below. Other interesting charge assignments with 
similar overall phenomenology are presented in Sec. IIVI 
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Table I: U(l)' charge assignments of the right-handed fermions and the scalar 4>. 



In the scalar sector, we adopt only one SM singlet scalar <f) with U(l)' charge 11. We can then write 
down the following renormalizablc couplings relevant for the neutrino masses 



C m = {mD)ijVL,iVR,j + \{ M R )ijl> Ri V Ri . 



S{vii.i + yi4> S^S 2 + yify S%S§ + S£S 7 + h.c, 



(3) 



where appropriate sums over i and j are understood. The rriD terms stem from electroweak symmetry 
breaking using the usual SM Higgs doublet H, while u>.; and yi are Yukawa couplings. Absorbing phases 
into the Sj we can take yj and one of the Wj to be real, while Mr can taken to be real and diagonal 
as well. Once (f) acquires a VEV, the ms matrix will be generated as mj = Wj(4>), leading to the MES 
structure in Eq. ([!} (discussed in detail in Sec. lIII A"|) . The other fermions 5*2,3,4.5,6.7 decouple from vl, vr 
and Si, and actually can be paired together to form three (stable) Dirac fermions ^1,2.3, to be discussed 
in Sec. ImCl 

Let us briefly comment on a theoretical constraint on the model. An inherent problem in any gauge 
theory involving abelian factors is the occurrence of a Landau pole, i.e., a scale at which the gauge 
coupling becomes so large that our perturbative calculations break down. In our model, the one-loop 
beta function of the U(l)' gauge coupling g 1 takes the form 



dln/i 



g' = 



16tt 2 



16tt 2 



(4) 



so the Landau pole of g' appears around the scale 



A, 



A' exp 



(5) 



where A' characterizes the U(l)' breaking scale. Inserting the U(l)' charges given in Tab. U we find 
b = 315, whereas for the 3 + 2 scheme from Tab. [XT] (which will be discussed later on in Sec. IIV A[) we 
have b = 75. For A' ~ 1 TcV and A L > AI Pl ~ 10 19 GeV, one obtains the constraints g'(A') < 0.08 
for the 3 + 1 case and g'(A') < 0.17 for 3 + 2 case. Alternatively, if we take the cutoff scale of the 
model to be the right-handed neutrino mass scale, i.e., A^ > Mr ~ 10 14 GeV, these bounds relax to 
g'(A') < 0.1 and g'(A') < 0.2 for 3 + 1 and 3 + 2, respectively. These upper bounds are stricter than the 
naive perturbativity bound g' 2 /4ir < 0(1)/ max(y') 2 . 



III. PHENOMENOLOGICAL CONSEQUENCES 



In this section we discuss the phenomenology of our new particles, with an emphasis on the novel 
effects in our framework. 



A. Neutrino Masses and Active— Sterile Mixing 



Let us consider the neutrino masses and active-sterile mixing in the U(l)' model. After the breaking 
of the gauge group SU(3) C x SU(2) L x U(1) Y x U{1)' to SU{3) C x C/(1)em, the full 13 x 13 mass matrix 
for the neutral fermions in the basis v = [y L 1; v L 2 , v L 3 , 1; v R 2 , v R 3 , S\, S%, S$, S%, , Sq, S?) reads 

(^MEs)rxr \ 
(Ms) 6x J 
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where the matrix Mmes reproduces the MES structure derived in Ref. @, Q (although our definition for 
ms differs by transposition), while Ms denotes the mass matrix of S2-7, explicitly given as 



Mi 



/ yi(<t>) \ 

yi(<p) 

y 2 (4>) 

y 2 (ct>) 

y 3 (<t>) . 

\ y a {4>) / 



(7) 



Obviously S2—7 decouple from the neutrino sector and can no longer be interpreted as right-handed 
neutrinos, because they do not partner-up with the SM neutrinos. The bare mass term Mr is unrestricted 
and can be large, as in the canonical seesaw case. We will consider this possibility here by setting 
Mr 3> mfl,ms, which leads to the effective low-energy neutrino mass matrix 

^4x4 _ _ f m D M^m% m D MR l ms\ 
" ~ ym^M- 1 mJ ) mgM^ms j ' 

for (vl, Sf)3 Such a mass matrix can be diagonalized by means of a unitary transformation as Mf, x4 = 
V diag(mi, m.2, m^, m 4 ) V T . Phenomenologically, the most interesting situation arises for to 5 3> mc, 
since the hierarchical structure of M 4xA allows us to apply the seesaw expansion once more, and arrive 
at the sterile neutrino mass 

to 4 ~ -rngM^ms , (9) 
together with the mass matrix for the three active neutrinos 

Ml x3 ~ -mDM^mJ) + moM^ms {mgM^ 1 m s )~ 1 mgM R 1 m^ ) = f7diag(mi, m 2 ,m 3 ) U T , (10) 
diagonalized by U. The 4x4 unitary mixing matrix V is approximately given by 

((1-^)17 R \ (n) 
\ -RW \-\b$r) 

with the active-sterile mixing vector R = tudM^ ms (rrioM^ ms) -1 = C(to_d/^s)- As a rough 
numerical estimation, for tud ~ 10 2 GeV, ms ~ 5 x 10 2 GeV and Mr ~ 2 x 10 14 GeV, one obtains 
the active neutrino mass scale m v ~ 0.05 eV, the sterile neutrino mass scale m s ~ 1.3 eV together with 
R ~ 0.2. This is in good agreement with the current global-fit data, i.e., |i?i| ~ 0.15 and Am^ ~ 1.8 eV 2 
for the 3 + 1 scheme [|| . 

For Yukawa couplings of order one, the observed large active-sterile mixing implies the scaling 
(4>)/{H) ~ 5-10. The new physics scale around TeV is hence not tuned to make LHC phenomenol- 
ogy most interesting, but comes directly from the neutrino sector. Actually — even though we obtain the 
magic TeV scale — the LHC implications of our model are rather boring, as we only expect small mixing 
effects in the Higgs and Z-boson interactions, to be discussed in the next section. 

Let us briefly comment on thermal leptogenesis [l3| in our framework. In principle, the additional 
singlet fermions may spoil the ordinary picture of leptogenesis since the right-handed neutrinos might 
predominately decay to sterile neutrinos instead of active neutrinos. This drawback can be easily cir- 
cumvented here by choosing the coupling of the lightest right-handed neutrino vr\ to the new states 
to be small, i.e., w\ <C 11*2,3. This will not modify the desired MES structure in the neutrino sector, 
but sufficiently increase the branching ratio of vr^i into SM particles, so standard thermal leptogenesis 
ensues. 



B. Bosonic Sector 

In this subsection we will briefly summarize the behavior of (j> and Z'. The scalar potential W in our 
model takes a simple form 

W = -i? H \H\ 2 + X H \Ht M 2 10| 2 + A^| 4 + 8 |ff| 2 |0| 2 , (12) 



§ On a more fundamental level, one can integrate out the heavy right-handed neutrinos vr at energies E -C Mr to 
generate the effective dimension-five Weinberg operators (mo)ij('mr))i e jLiHH^ L^/((H) 2 (Mn)jj), <j> 2 SiS1/(Mr)h 
and (mjj)ijWj LiHSi<j>^ /((H) (Mn)jj), which were the starting point in Refs. [id . 
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where H denotes the SM Higgs doublet, and can be decomposed as 4> = (Re <j) + i Im <f>) / v2 = ((Re0) + 
f + i Im </>) / y/2. After symmetry breaking and in unitary gauge, Im is absorbed by the Z' boson, giving 
it a mass Mz> = \ Wg'{<j>)\. Due to the 5 term in the scalar potential, we have a generic mixing between 
the remaining real field ip and the neutral scalar h contained in H , i.e., 



cos 9 — sm 
sin 9 cos 9 




(13) 



where h\ and h 2 are the physical mass eigenstatcs, and the mixing angle 9 is given by 



sin 29 = MB . (14) 



A non-zero 5 — and hence — opens the well-known Higgs portal [14j for the DM production/annihilation, 
which will be discussed in the next section. 

The Higgs portal |</>| 2 |-ff| 2 aside, there is one more renormalizablc, gauge-invariant operator that will 
induce a coupling between the SM and DM sectors, namely the kinetic-mixing operator sin£ Fy ' F'^ [l5jj. 
This off-diagonal kinetic term involving the hypercharge and U(l)' field strength tensors will induce a 
coupling of the physical Z' boson to the hypercharge current. The relevant phenomenology of the resulting 
interaction between the SM and DM particles can be found in e.g. Refs. (m Il7j. 



C. Dark Matter 



As we mentioned before, the singlet fermions 5*2,3,4,5,6,7 are the DM candidates in our model. To see 
this point, we write down the full Lagrangian for the right-handed singlets 5*2 and S3: 



£s 2 , 3 = iS 2 y<(d, - i(-5g')Z'JS 2 + iS 37 M (d M - i(-6g')Z'JS 3 + yi(<t>T 3 S 2 + h.c.) 



(15) 



By defining the Dirac field vPi = S 2 + S3, the above Lagrangian can be rewritten in unitary gauge as 



C S2 3 =1*17^*! +g'Z'V 1 ^ 



(-5)- (-6) , (-5) + (-6). 



2/1 



75 ) *i + -J|((Re(^))+vj) 



(16) 



After spontaneous symmetry breaking, \&i acquires a Dirac mass Mi = -yi(Re<f>)/V2. Similarly, we can 
define ^2 = S4 + S5 and ^3 = Sq + Sy for S4,5 and Sg,7, and obtain altogether the DM Lagrangian 



C 



DM 



= E 

J=l,2,3 



Ma 



(17) 



where we defined Yj = Y'(Sj). The stability of these fields will be discussed below, but let us first take 
a look at the interactions involving the will-be sterile neutrino Si, given by the Lagrangian 



C Sl = iS l7 ^ - i(ll 5 ')^)Si + M f SfPR,i + h.c.) 



(18) 



The important part is the Z' interaction, as it allows for the annihilation — > Z' — > S1S1. Since the 
physical sterile neutrino v s = consists mainly of Si, but contains a not-too-small part of the active 
neutrinos v e ^ r , this process connects the DM to the SM sector. Specifically, this "neutrino portal" takes 
the form 



C 



^ — portal 



~2 P 



- ll75)*i + ^27^(13 - H75)*2 + *37 M (-7 + ll75)*3 



11 KVij tpirfitvi +^1%) 

«,i=i 



(19) 



where the four light mass eigenstates Vj are written as Majorana spinors and the unitary matrix V is 
defined in Eq. (|TTj) . 
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Figure 1: Visualization of the different scales in our framework, as well as the relevant interaction channels. 
The red connection between the DM and SM sectors represents the well-known kinetic-mixing and Higgs portals 
(based on vector and scalar mixing respectively), while the blue interactions are relevant for the neutrino portal 
(based on fermion mixing). Interactions with the vrj are highly suppressed and not shown. 



Wc further note that the heavier DM particles can also convert to the lighter ones, i.e., "Fi'Fi —> 'Fj'Fj 
via the s-channel exchange of Z' or 4>. Moreover, ^ may also annihilate to Z' and </>, which can enhance 
the total annihilation cross section significantly. 

The model content and relevant scales are illustrated in Fig. [TJ The (self-interacting) DM sector couples 
to the SM just like all models with a dark symmetry U (1)dm, namely through scalar mixing (Higgs portal, 
parameterized through 5) and vector mixing (kinetic- mixing portal, parameterized through £). However, 
due to the gauge interactions of the DM with the sterile neutrinos, a new portal through fermion mixing 
(neutrino portal) opens up in our model. Since this portal is not often discussed in the literature (see 
however Refs. [Ig,[l9]), we wm " f° cus 011 it m the remainder of this paper. 



1. Stability 



It is obvious that the $ \ fields in Eq. (fTTj) are stable, since there exists an accidental global U(l) 3 
symmetry shifting the phases of 'F,. The occurrence of several stable DM particles (multicomponent 
DM) results in numerous interesting effects (see Refs. [20j | for some early work). The underlying reason 
for the stability in our case is the remaining exact Zn symmetry after the spontaneous breakdown of the 
U(l)'. The \&j form representations under this discrete gauge group with charges 6, 1 and 2 (modulo 11), 
which stabilizes at least the lightest of them, even when higher-dimensional operators are considered. 

While our model is renormalizable, we expect it to be only valid up to a certain cutoff scale A, either 
because quantum gravity takes over, or because sooner or later we will hit the U(\)' Landau pole — 
as discussed in Sec. |TTJ At the cutoff scale, higher-dimensional operators might be generated, and in 
our models these will always include <p 2 SiSf/A and the Weinberg operator [2l[ for i/£-Majorana masses. 
Taking A ~ Mpi does not destroy the discussed MES structure if (4>) < 10 TeV. For the charge assignment 
here, there are also dimension-six operators like 64 fig S^/ A 2 , which break the global U(l) 3 to a U(l) 
symmetry, so only one stable Dirac fermion survives. However, since these operators are highly suppressed 
for A ~ Mpi, the resulting lifetimes are typically longer than the age of the universe, and thus we will 
not include them in our discussions below, but take all three "F., to be independently stable. 



2. Relic Density and Thermal History 



We will now discuss the interplay of the three portals (Higgs-, kinetic-mixing-, and neutrino portal) and 
identify some valid regions in the parameter space where the correct relic density for "F^ can be obtained. 
Note that the mixing parameters 5 and £ are the only new physics parameters we assume to be small in 
this paper, all other couplings are somewhat "natural" . We restrict ourselves to small mixing parameters 
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solely for simplicity, as larger values lead to very constrained effects, see Refs. [HI, [H, [13, H3- A more 
detailed parameter scan (including thermal leptogenesis) will be performed in a separate publication [23| . 

We only consider freeze-out scenarios. Also note that we always end up with a thermalized sterile 
neutrino at the epoch of neutrino decoupling, so the usual cosmological bounds on N e g and ^2 vn v hold . 
This is to be expected in models with light sterile neutrinos, and can be solved on the astrophysics side, 
or by choosing a smaller-than-eV mass for the sterile neutrino. 

Case A: 6, £ = 0. To check the validity of the neutrino portal, we first turn off the Higgs- and kinetic- 
mixing portals by setting 6 = £ = (or at least small enough to be negligible). In this case, the only 
connection between the new physics sector and the SM comes from active-sterile mixing, or, at a more 
fundamental level, from the exchange of heavy right-handed neutrinos. Integrating out the vr gives for 
example the operator LHSi(f>/Mn (using order one Yukawa couplings), which gives a rough scattering 
rate for LH 4+ Si</> around ~ T 3 /Mj| — to be compared to the expansion rate in the early universe 
~ y/g^T 2 /Mpi — which puts all particles in equilibrium above T > 10 10 GeV.^ Below that temperature, 
the two sectors SM and DM (consisting of Z', <j> and Sj) evolve independently, while the temperature 
decreases due to expansion of the Universe in both sectors. Nothing really happens until T ~ TeV, when 
the tyj freeze-out occurs. To deplete the tyj abundance fast enough, we can make use of the neutrino 
portal, i.e., the annihilation of the lightest tyj into v s around the Z' resonance (assuming that the heavier 
tyj annihilate sufficiently fast into the lightest "J^-, see Fig. [5] for illustrations). 

After freeze-out, we have overall three decoupled sectors SM, tyj and v s — all with different tempera- 
tures. Above (active) neutrino decoupling, the Universe was radiation dominated, so only the temperature 
of v s and the relativistic degrees of freedom in the SM sector are of interest and will be calculated now. 
Using conservation of entropy in the two sectors SM and DM, we have the equalities 



7 f M T 3 a 3 



= 9 



SM 



and 



DM T 3 a 3 



DM T 3 
9* 1 DM' 



(20) 



where denotes the effective number of relativistic degrees of freedom in sector X, a the scale factor, t scp 
the time when the two sectors just separated from equilibrium (i.e., at temperatures around 10 10 GeV) 
and tf the final time we are interested in, namely close to active neutrino decoupling (e.g. when Tqm ~ 
lOMeV). At tf, the SM sector consists of photons, electrons and neutrinos, while the DM sector only 
has the relativistic Si ~ v s , so wc find 



SM 



g. UM (t scp ) gt M {tf) 
9? M (tf) 9l M (t scp ) 



1/3 



/ 65/4 43/4 \ 1/3 
V 7/4 427/4 J 



0.98. 



(21) 



Ignoring active-sterile oscillations, this would make the sterile neutrinos slightly colder than the active 
ones at decoupling, alleviating cosmological constraints to some degree (the one sterile neutrino effectively 
contributes only AN e g = (T„,/Tsm) 4 — 0.92 additional neutrinos to the energy density). However, for the 
sterile neutrino parameters relevant for the short-baseline anomalies, i.e. m s ~ eV, S ~ 0.1, active-sterile 
oscillations will become effective around T ~ 100 MeV-1 MeV |24| , once again connecting the SM bath 
and v s and thus thcrmalizing the sterile neutrino at neutrino decoupling. Note that the usual discussions 
of active-sterile oscillations at these temperatures are not readily applicable, as our model starts with 
abundant v s and self-interactions mediated by Z' (freezing out around Tdm ~ 10 MeV). In any case, the 
cosmological bound on relativistic degrees of freedom is expected to be approximately valid in our model 
and will be discussed in more detail in a separate paper [23| . 

Case B: 5^0. Let us open the Higgs portal. The thermal evolution is similar to case A, but values 
S > 10 -7 will put <j> in equilibrium with the SM at temperatures below T ~ 10 TeV, because the scattering 
rate hh <h> <p(f> goes with 5 2 T/Att [l?]]. 4> an d the rest of the DM sector (Z 1 , "J j and v s ) are in equilibrium 
through U(l)' gauge interactions (for not too small gauge coupling g'), so SM and DM are in equilibrium 
around DM freeze-out. For the freeze-out we can again use the neutrino portal, i.e., resonant annihilation 
, I n I' — > Z' — > v s v s . As (f> and Z' go out of equilibrium around the same time, the connection between 
the SM sector and v s is severed and the two evolve independently for a while, until they are reconnected 
around T ~ 10 MeV by active-sterile neutrino oscillations. 

In a different region of paramctcrspacc, we can make use of the resonant annihilation of DM into SM 
particles via scalars, i.e., the Higgs portal in the way it is intended. The discussion is then completely 
analogous to other C/(1)dm models, so we refer the interested reader to Ref. [22j for a recent evaluation. 



^ We assume a sufficiently high reheating temperature after inflation. 
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Figure 2: fi* /i 2 versus the Z' mass M z / for degenerate (top panels) and hierarchical (middle and bottom panels) 
DM masses. The VEV was set to (</>) = 1.5 TeV and the scalar mass is labeled on the plot. The red, green 
and blue lines show the relic density of ^i, ^2 and $ 3 respectively, while the black dashed line gives the full 
fi*/i 2 = Qigjh 2 . The horizontal pink band represents the relic density measured by WMAP [25| (la range). 



Case C: £ 7^ 0. A very similar discussion can be made for an open kinetic-mixing portal. Again small 
values £ > 10~ 7 suffice to reach thermal equilibrium of the SM and DM sectors, e.g. through scattering 
Zh <H- Z'h. The thermal evolution then closely resembles that of case B, with some minor differences: 
The Z'-Z mixing couples v s to the SM, so Z' interactions keep v s thermalized a while longer before it 
decouples and finally reconnects with the SM. Furthermore, the DM annihilation around the Z' resonance 
contains a small branching ratio into SM particles. 

In order to illustrate the feasibility of the DM candidates via the neutrino portal, we implement the 
model in micrOMEGAs [26| and evaluate the relic density of DM particles ^i. Here we assume 5 and 
£ to be large enough to thermalize the SM and DM sectors at DM freeze-out, but small enough to be 
negligible in the numerical calculation of the relic density, i.e., we only consider DM freeze-out using the 
neutrino portal. The scalar VEV is taken to be ((f) = 1.5 TeV as an example. The gauge coupling g' 
are therefore obtained from the relation of Z' mass and ((f)) . As shown in the upper panels of Fig. [2j 
a resonance appears at ~ Mz>/2, and the relic density f2^/i 2 ~ 0.1 measured by WMAP [25| can 
be achieved. In the degenerate case (i.e., Mi ~ M2 ~ M3), the 'Fi contribution to is dominating 

because it has the smallest Z' coupling. Moreover, in case of a small scalar mass, e.g., = 500 GeV, a 
new channel ^"F — > Z'(f> is open for light Z', which is observed from the upper-right panel of Fig. [2] For 
the case of non-degenerate spectrum (i.e., M\ 7^ M-x 7^ M3), the most significant contribution to the relic 
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Figure 3: Coupling of DM to leptons. 



density may come from either ^i, ^2 or ^3, depending on the specific fermion spectrum as well as the 
scalar and vector masses. As can be seen in the middle and lower panels of Fig. [3J the >Fi contribution 
to the relic density typically dominates, but there exist model parameters that make ^2 or ^3 the main 
DM particle. 

3. Direct and Indirect Detection 

The neutrino portal discussed so far does not lead to any direct detection signals, because the cross 
sections are highly suppressed. Loop processes connecting >F to SM fermions, e.g. as in Fig. [3J van- 
ish in case of degenerate active-sterile masses, so these amplitudes are suppressed by tiny factors like 
ATO^/O(100GeV) 2 ~ 1CT 22 . 

Indirect detection might naively be more fruitfull, because the annihilation of the 'Fj in the Galactic 
Center or halo leads to two back-to-back neutrinos with energies ~ Mj (whichever >Fj is sufficiently 
abundant), which is an ideal signal for neutrino telescopes like IceCube.** However, since we considered 
<Fj to be a thermal relic, the self-annihilation cross section is set by the relic density, which is too small 
to be probed [I?} — even though the branching ratio into neutrinos is ~ 100%, so the signal is as clear as 
it gets. 

Direct and indirect detection measurements are of course sensitive to the Higgs- and kinetic-mixing 
portal parameters S and £, as discussed in the literature (e.g. Ref. [l9[ discusses the Higgs portal in a 
framework similar to the neutrino portal). 

IV. OTHER INTERESTING CHARGE ASSIGNMENTS 

Having focussed on one specific example using the charges from Tab. U we will now briefly present 
other charge assignments with interesting phenomenology. In all cases, we only introduce one additional 
scalar </>, so the results concerning scalar and vector interactions remain unchanged — different numerical 
values for the charges aside. Only the sterile neutrino and dark matter sector will be slightly modified. 

A. More Light Sterile Neutrinos 

The introduction of n > 2 light sterile neutrinos (3 + n scheme) increases the number of new parameters 
and most importantly allows for CP-violation in the effective oscillation analysis (28|. This feature can 
significantly improve the fit to neutrino oscillation data and has been studied extensively 0, 0, [H[ . Note 
that the tension with standard ACDM cosmology typically worsens, depending on the used data sets [30j . 

We can easily modify the above U(l)' framework to accommodate the 3 + 2 MES scheme, by choosing 
different charges for the ten right-handed neutrinos: we need at least one more neutral vr^ to generate 
the necessary light mass squared differences. Now we have to find charges that treat two of Si the same 
(without loss of generality Si and S2), i.e., Y'(Si) = F'^), so these will become our two light neutrinos 
after coupling them to a scalar (f>. We can once again find exotic charges in such a way that the decoupled 
Sj become massive by coupling to the same scalar, the magic number for this to happen seems to be six. 
See Tab.|TT]for a valid anomaly-free charge assignment with the desired properties (also used in Ref. [lol]). 
After breaking the U (1)' and the electroweak symmetry, the 13 x 13 mass matrix for the neutral fermions 



* The DM-micleon cross section in our model is too small to efficiently capture DM inside the Sun or Earth, so we have to 
rely on astrophysical objects with high DM density. 



10 





VR,l 


VR,2 


VR,3 


VR,4 


Si 


s 2 


s 3 


5 4 


s 5 


Se 





Y' 














-5 


-5 


-1 


6 


2 


3 


5 



Table II: Exotic U(l)' charge assignments of the right-handed fermions and the scalar <f> to obtain the 3 + 2 MES 
scheme. 



takes the form 

m= ^Mmes) 9x9 i (22) 

where the 9x9 matrix .Mmes m the basis (y e , v^, v T , v c R 1; z/j^ 2 , z/^ 3 , z/^ 4 , SI, S%) is the obvious extension 
of the MES structure from Eq. ([!) for the 3 + 2 scheme, while Ais denotes the mass matrix of S3-Q, i.e., 




Mi 



( Vl {4>) \ 

vi(4>) 000 

y 2 {4>) 

V ite<0> y 



(23) 



resulting in two Dirac fermions, decoupled from the neutrino sector. 

Compared to the 3+1 scheme discussed so far, the scalar sector is identical, whereas the dark matter 
sector is slightly modified because we have only two stable Dirac fermions (protected by the remaining 
discrete gauge group Z5) instead of three, but two light sterile neutrinos instead of one. This does not 
influence the qualitative behavior significantly. 

The expressions from Sec. IIII Al for the neutrino masses go through in the same manner, i.e., we still 
have 

M 3x3 _ -m D M- 1 m^ ) + moM^ms {mlM- 1 m s )~ 1 m T s M- 1 m T D , M 2 V * 2 ~ -m^M^ms , (24) 

for the masses, where we assumed mo <C rris <C Mr, and run, m$ and Mr are 3 x n, n x 2, and n x n 
matrices, respectively. The number n of gauge group singlets i/r should be at least 4 to make enough 
light neutrinos massive, n > 5 makes all five light neutrinos massive. The active-sterile mixing is again 
0(mD/iris), so the required values O(0. 1) put the U(l)' breaking scale naturally in the TeV range. 

Let us briefly comment on the thermal evolution of the universe in this model. Seeing as the number 
of degrees of freedom is smaller (larger) at t scp (tf) compared to the 3+1 scheme of Sec. IIII C 21 the 
sterile neutrino bath is colder than the SM bath (prior to neutrino decoupling) by a factor of ~ 0.75. 
Without active-sterile neutrino oscillations, this would mean that the two sterile neutrinos effectively 
only contribute AN e {{ ~ 0.6 additional neutrino species to the energy density, alleviating cosmological 
bounds. It is of course to be expected that active-sterile oscillations before neutrino decoupling generate 
thermal equilibrium among the neutrinos, giving rise to the usual constraints. 

For completeness, we give an assignment for the 3 + 3 case, which has been fitted to the neutrino 
anomalies in Ref. [3l|. To make at least five light neutrinos massive, we need five Vr. A possible charge 
assignment for nine Sj is then (7, 7, 7, 2, —9, —1, —6, —4, —3), with one scalar 0^7. This leads to three 
light sterile neutrinos and three stable Dirac DM particles (protected by the remaining Z7). 



B. Majorana Dark Matter 



Having focused on Dirac DM in the main text for no particular reason, we will now give an example 
with Majorana DM. For the 3 + 1 MES scheme, we take the exotic charges (6, —3, —3, 2, —8, —1, 7) for 
the Si and one scalar with charge Y'((j>) = 6. The VEV of (j> breaks U(l)' — > Xq, S\ will again become 
the sterile neutrino, while S2 and S3 share the most general Majorana mass matrix — which we can take 
to be diagonal without loss of generality — resulting in two Majorana fermions ^1,2- (S4, S§, Sq, S7) share 
the mass matrix 



Mi 



( 


















V 

resulting in two Dirac fermions ^3,4; all 
representations ^1,2 ~ 3 ~ (1,0), ^3 ~ 2 ~ 
the mass spectrum, we can obtain a stable Majorana fermion. 






2/2(0) 




(25) 



are decoupled from the neutrino sector. These particles form 
(0,2), ^4 ~ 1 ~ (1, 1) under 7L§ = Z2 x Z3, so depending on 
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C. Unstable Dark Matter 

The charges for the Si and <p discussed so far have been chosen in such a way that the spontaneous 
breaking of U(l)' leaves a nontrivial Z^v that stabilizes the DM candidates. This is of course not a 
generic feature of exotic charges, but just a convenient choice to obtain exactly stable particles. Let us 
briefly comment on unstable DM candidates: Taking (1, —10, 9, —7, 6, —11, 12) for the 3 + 1 scheme with 
a scalar <f> ~ 1 gives three Dirac DM candidates — S2 + S$, $>2 = S4 + S£, ^3 = Sq + Sj, which are 
independently stable due to an accidental global t/(l) 3 symmetry. However, with this charge assignment, 
there is no leftover Zjv symmetry protecting this stability. Similar to the discussion in Sec. MI CTl we can 
discuss higher-dimensional operators. For the charge assignment here, there are already dimension-five 
operators 

4> 2 S 3 S C 4 /A, <f 2 S C 3 S 6 /A, 4> 2 S 2 S C 7 /A, (26) 

which break the global t/(l) 3 to a U(l) symmetry, so only one stable Dirac fermion survives. Even this 
stability is not exact, as there are operators like 4> & S b vn/A 5 which break the global U(l) and lead to 
DM decay. In this particular example — and for A ~ Mp\ — the decay would be suppressed enough to 
still allow valid DM, but in principle there are charge assignments with decaying DM, or even no DM 
candidate at all. 



V. CONCLUSION 

Generating small sterile neutrino masses via the same seesaw mechanism that suppresses active neutrino 
masses requires a specific structure in the neutral fermion mass matrix. We showed how this so-called 
MES structure can be obtained from a new spontaneously broken U(l)' symmetry, under which the 
"sterile" neutrino is charged. Heavily mixed eV-scale steriles hint at a U(l)' breaking scale around TeV. 
Additional anomaly-cancelling fermions need to carry exotic U(l)' charges in order to not spoil the MES 
structure, which coincidentally stabilizes one or more of them (without the need for additional discrete 
symmetries). The main connection between this multicomponent dark matter sector and the Standard 
Model is the active-sterile mixing (neutrino portal). We discussed how the dark matter annihilation 
almost exclusively into sterile neutrinos can be used to obtain the measured relic density, and also the 
interplay with the other two portals (Higgs- and kinetic-mixing portals). 

We focussed on a few specific examples, but the presented framework of exotic charges obviously 
provides a rich playground for model building, depending on the used charges and number of new particles. 
Worthwhile extensions with U(l)' charged SM fermions, e.g. B — L type symmetries, can be obtained 
with slightly more complicated scalar sectors and will be discussed elsewhere.' ' 
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